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ABSTRACT 

Variants of mitochondrial DNA (mtDNA) are 
commonly used as markers to track human evolu- 
tion because of the high sequence divergence and 
exclusive maternal inheritance. It is assumed that 
the inheritance is clonal, i.e. that mtDNA is 
transmitted between generations without germline 
recombination. In contrast to this assumption, a 
number of studies have reported the presence of 
recombinant mtDNA molecules in cell lines and 
animal tissues, including humans. If germline 
recombination of mtDNA is frequent, it would 
strongly impact phylogenetic and population 
studies by altering estimates of coalescent time 
and branch lengths in phylogenetic trees. 
Unfortunately, this whole area is controversial and 
the experimental approaches have been widely 
criticized as they often depend on polymerase 
chain reaction (PCR) amplification of mtDNA and/ 
or involve studies of transformed cell lines. In this 
study, we used an in vivo mouse model that has had 
germline heteroplasmy for a defined set of mtDNA 
mutations for more than 50 generations. To assess 
recombination, we adapted and validated a method 
based on cloning of single mtDNA molecules in the 
X phage, without prior PCR amplification, followed 
by subsequent mutation analysis. We screened 2922 
mtDNA molecules and found no germline recombin- 
ation after transmission of mtDNA under genetically 
and evolutionary relevant conditions in mammals. 

INTRODUCTION 

Mutations in mitochondrial DNA (mtDNA) are the most 
frequently used genetic marker in phylogenetic and 



population studies in mammals. This popularity rests 
upon three pillars: the high sequence divergence of 
mtDNA, the exclusive maternal transmission and the 
assumed lack of germline recombination. This means 
that mtDNA is used as a molecular marker showing 
patterns of clonal inheritance, which reduces the complex- 
ity of the analyses as bi-parental inheritance does not have 
to be taken into account. The possible occurrence of 
mtDNA recombination in the mammalian maternal 
germline has far reaching implications, because even 
modest frequencies of germline recombination would sub- 
stantially affect coalescent time estimates and branch 
lengths in phylogenetic trees (1). During the last two 
decades numerous reports have questioned the assump- 
tion that there is no mtDNA recombination based on 
analyses of human patient tissues (2,3), tissue culture 
cells (4,5) and tissues from interspecies animal hybrids 
(6,7). However, many of these studies have relied heavily 
on polymerase chain reaction (PCR) amplification with no 
or little control for recombination by jumping PCR, which 
may produce chimeric PCR products (8) that falsely are 
interpreted as evidence for in vivo recombination. Another 
common limitation of these reports is that they mostly 
cover events in somatic tissues or are based on analyses 
of transformed cell lines and therefore have little relevance 
for understanding mtDNA evolution of animals. To com- 
plicate matters further, it should be noted that bi-parental 
inheritance of mtDNA has been robustly documented in 
some organisms, like the mussel Mytilus galloprovincialis 
(9). Within mammals, heteroplasmy, i.e. the presence of 
two mtDNA variants within the same individual, is more 
common than previously thought (10,11) and paternal 
leakage of mtDNA has been reported under unusual 
conditions, such as interspecific mouse crosses (12,13). 
Hence, there are certainly documented circumstances 
where different types of mtDNA molecules coexist in an 
individual, but it is unclear to what extent recombination 
of mtDNA actually occurs in the mammalian germline. 
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Unfortunately, technical limitations combined with the 
lack of appropriate in vivo experimental systems has left 
the research field in uncertainty, with claims both for and 
against recombination of mammalian mtDNA being ques- 
tioned. In this report, we have used a maternal mouse line 
that has been heteroplasmic for a defined set of mtDNA 
mutations for more than 50 generations in combination 
with a protocol based on X phage cloning for sequence 
analysis of single mtDNA molecules. We analyzed 2922 
mtDNA molecules and found no recombinant molecules 
when mtDNA was transmitted in the maternal germline 
under genetically and evolutionarily relevant conditions. 
These findings thus demonstrate that recombination does 
not contribute to mtDNA variation in mammals. 

MATERIALS AND METHODS 

Mouse breeding 

The NZB-BALB/c heteroplasmic mouse line was 
generated by cytoplast transfer, as previously described 
(14), and has been maintained by backcrossing of 
heteroplasmic females to inbred BALB/c males for more 
than 50 generations. 

Preparation of mtDNA for X phage cloning 

Liver and kidney mtDNA was extracted and purified 
according to standard protocols. In brief, tissues were 
homogenized with 15 strokes of a teflon homogenizer at 
1500rpm in mitochondrial isolation buffer [MIB; 225 mM 
mannitol, 7.5 mM sucrose, 1 mM 4-(2-hydroxyethyl)-l- 
piperazineethanesulfonic acid (HEPES), 1 mM ethylene- 
diaminetetraacetic acid (EDTA) pH 8.0, 0.1% bovine 
serum albumin and 0.025% (3-mercaptoethanol]. Cell 
debris was pelleted by low speed centrifugation (lOmin 
at 600 g) and mitochondria were recovered from the super- 
natant by high-speed centrifugation (lOmin at 5000g). 
Mitochondrial pellets were washed three times in MIB, 
using a loose-fit glass homogenizer to resuspend the 
pellets (Wheaton, USA). The final pellet was placed on a 
two-step sucrose gradient cushion (1.5M/1M), followed 
by lh centrifugation at 77 000 g (SW32Ti, Sorvall). 
Mitochondria were recovered from the gradient, pelleted 
in mitochondrial gradient buffer (MGB; 10 mM 
HEPES pH 7.8, lOmM EDTA pH8.0), re-suspended in 
mitochondrial lysis buffer (MLB; 75mM NaCl, 50 mM 
EDTA pH 8.0, 20 mM HEPES pH 7.8) and lysed with 
1% Sarkosyl, followed by phenol/chloroform extraction 
and isopropanol/NaCl precipitation. The concentration 
mtDNA was determined on a DyNAQuant 200 (Hoefer). 

X Phage cloning 

Sucrose-gradient purified mtDNA (100-200 ng) was 
digested with the Bglll endonuclease (NewEngland 
BioLabs). Appropriate cohesive ends were generated by 
Klenow DNA polymerase extension in the presence of 
0.17mM dGTP and dATP, followed by standard 
phenol/chloroform extraction and ethanol/salt precipita- 
tion supplemented with 1 ul Pellet paint (Novagen). 
The DNA pellet was dissolved in 1.8 ul ddH 2 0, 1 ul 1 A. 
phage vector arms (Stratagene), 0.35 ul lOx T4 DNA 



ligase buffer, 0.3 ul T4 DNA ligase (NewEngland 
BioLabs) — with a DNA concentration of >0.2ug/ul — 
and incubated for 20 h at +4°C. The ligation reaction 
was added to an in vitro packaging extract (Stratagene) 
and incubated for 90 min at room temperature. 
Packaged phage library titers were prepared in 500 ul of 
SM buffer (lOmM NaCl, 8mM MgS0 4 , 50mM TRIS 
HC1 pH 7.5) and 20 ul of chloroform. Individual X clones 
were expanded in an infected host Escherichia coli strain, 
which had entered the lytic cycle for plaque formation. 
The host E. coli was grown according to conditions in 
the provided protocol and pelleted and dissolved in 
10 mM MgS0 4 to an OD 600 of 0.5. Two hundred ul of 
E. coli cells were then incubated with an appropriate 
volume of phage (0.01-2 ul) at +37°C for 15 min. The 
phage — E. coli solution was then diluted with 3.5 ml of 
liquid NZYM top agar, spread on NZYM agar plates 
and incubated over night at +37°C. The following day 
individual plaques were picked using a 1-ml pipette tip 
and eluted into a suitable solution for further analysis. 
Re-plating experiments suggest that a single plaque 
contains approximately 50000 copies of the cloned 
mtDNA. 

PCR and restriction fragment length polymorphism 
analysis 

One microlitre of eluted phage was used as PCR template 
per reaction. DNA was extracted from phage particles by 
a 10-min incubation at 95° C as the initial step in the PCR. 
PCR primers were designed to amplify appropriate parts 
of the mtDNA genome to be used in the restriction 
fragment length polymorphism (RFLP) screening for 
recombination (Supplementary Table SI). PCR products 
were cleaved with the appropriate endonucleases 
(Supplementary Table S2) and the digested PCR 
products separated by gel electrophoresis in 1.5% 
agarose gels. 

PCR and sequencing of full-length mtDNA molecules 

PCR primers, covering the entire mtDNA sequence 
(Supplementary Table S3) and containing M 13 sequencing 
tags at their 5'-ends, were designed to amplify 27 
overlapping PCR products. PCR products were purified 
with AMPure (Agencourt) and sequenced using BigDye, 
version 3.1 chemistry (Applied Biosystems). Cycle 
sequencing reactions were purified using X-terminator 
solution (Applied Biosystems) and run on a 3130xL capil- 
lary sequencer. Sequences were aligned to the C57BL/6J 
reference sequence (accession NC 005089.1) using the 
Seqscape V2.5 software (Applied Biosystems) and all poly- 
morphisms were numbered according to the C57BL/6J 
reference sequence. 

Quantification of heteroplasmy levels 

Relative heteroplasmy levels for whole-tissue extracts, 
isolated mitochondria and phage titer were quantified by 
last-cycle fluorescent PCR-RFLP analysis. The relative 
heteroplasmy levels of the final clones were also quantified 
by PCR-RFLP and the cleaved products were separated 
on 1.5 % agarose gels. PCR and RFLP analyses were 
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performed as previously described (15), using primers con- 
taining a 5' M13-tag and an internal digestion control 
(forward: TGT AAA ACG ACG GCC AGT TCT CCG 
TGC TAC CTA AAC ACC; reverse: CAG GAA ACA 
GCT ATG ACC TAG GTT GGT GCC GGA TAT TGT 
G). For the last-cycle fluorescent analysis, the final cycle 
of a standard PCR was spiked with 5uM M13F-hex 
primer. PCR products were purified using Agencourt 
Ampure (Beckman Coulter) and an aliquot was digested 
by NlalV endonuclease (New England Biolabs). The 
digested product (1 ul) was separated by fragment 
analysis on a 3130xL DNA Analyzer and analyzed using 
the GeneMapper v. 4.0 software (Applied Biosystems). 
Peak heights were used to determine relative heteroplasmy 
levels. Average values from three independent experiments 
were used for the quantification of mtDNA mutation 
levels. 

RESULTS 

PCR amplification induces abundant in vitro 
recombination artifacts 

The amplification of DNA by PCR can lead to the 
creation of artificial chimeric molecules by polymerase 
jumping (8). We have previously applied a standard 
PCR-clone-sequence protocol to detect mtDNA muta- 
tions in C57BL/6N mice (16,17). From these available 
datasets, we extracted sequences that had been filtered 
due to similarity with nuclear mtDNA pseudogenes 
(NuMts). We identified NuMts in 41 of 3940 clones 
(~1% of samples). Strikingly, nine of the clones 
(Figure 1) showed evidence of eight independent PCR- 
induced recombination events between NuMts and 
mtDNA (Figure 1). This high frequency was present 
despite the very low abundance of nuclear DNA 
relative to mtDNA in the samples. Cloning and subse- 
quent sequencing of PCR products is therefore not 
suitable for identifying mtDNA recombination events 
because artificial in vitro generated chimeric molecules 



are very frequent and indistinguishable from true in vivo 
recombinants. 

Cloning of mtDNA in the X phage does not create 
artificially chimeric mtDNA molecules 

The unsuitability of a PCR-clone-sequencing approach to 
study mtDNA recombination prompted us to develop an 
alternative method to capture large stretches of native 
mtDNA. The X phage has a long tradition of being the 
vehicle of choice to clone a variety of genomic targets. 
Utilizing this vector, we developed a direct cloning 
protocol that enables us to trap single full-length 
mtDNA molecules from mitochondrial-enriched DNA 
preparations using as little as 10 ng of input DNA. By 
sequencing entire phage inserts we were able to show ex- 
perimentally that X phage cloning can successfully trap 
full-length mtDNA molecules. Next, we created an artifi- 
cially heteroplasmic sample by mixing liver tissue from 
C57BL/6 and NZB mice, that differ by 90 positions in 
their mtDNA, to experimentally test whether chimeric 
molecules are artificially created by the X cloning 
approach. A standard mtDNA extraction procedure was 
then performed on the mixture of liver tissue from the two 
mouse strains, followed by X phage cloning. In order to 
screen a large number of mtDNA clones, we chose 
four strain-specific polymorphisms distributed across the 
mitochondrial genome (nucleotide positions 2814, 5463, 
9985 and 15 657) for RFLP analysis (Figure 2a and b). 
Individual plaques were analyzed for recombinants 
carrying both NZB and C57BL/6 mtDNA variants. In 
total, we screened 100 molecules and detected no recom- 
binant molecules. 

Cloning of mtDNA in the X phage is unbiased 

In order to ensure that the cloning method is not biased 
toward a subgroup of mtDNA molecules, we extracted 
mtDNA from a previously published mouse line with a 
heteroplasmic 3875delC mutation in the tRNA Met gene 
(15). We scored the 3875delC genotype in 830 cloned 
full-length mtDNA molecules from two heteroplasmic 
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Figure 1. Sequence of artificially created recombinant mtDNA molecules. The sequence patterns of mtDNA molecules analyzed by the PCR-clone- 
sequencing method are shown. The mtDNA reference sequence variants are depicted in white and the NuMts in grey. In total, 9 out of 41 
of clones showed evidence of PCR-induced recombination between NuMts and mtDNA. Asterisk indicates two clones generated by the same 
in vitro event. 
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Figure 2. RFLP analyses of mtDNA in heteroplasmic mice, (a) RFLP patterns obtained from restriction enzyme digestion of X phage cloned 
mtDNA molecules from liver showing pure BALB/c (lanes 1^1) or NZB (lanes 9-12) patterns. One cloned molecule (lanes 5-8) shows a mixed 
NZB/BALB/c cleavage pattern. The screened nucleotide variants (2814, 5463, 9985 and 15 657) are shown on top of each lane, (b) Schematic 
depiction of mouse mtDNA with the non-coding control region (yellow) and genes encoding proteins (green), ribosomal RNAs (rRNAs; blue) and 
transfer RNAs (tRNAs; brown). The four sites that were used for screening to distinguish between NZB and BALB/c mtDNA (nucleotide variants 
2814, 5463, 9985 and 15 657) are indicated with arrows. 



Table 1. The A phage protocol does not favor cloning of specific 
mtDNA genotypes 
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tRNA Met mice (418 and 412 for mice 1 and 2, respect- 
ively), by PCR-RFLP. Heteroplasmy levels in liver 
tissue, mtDNA extracts and the phage titers for the 
two respective animals were quantified by last-cycle fluor- 
escent PCR-RFLP analysis. Subsequent comparison of 
heteroplasmy levels in the liver tissue, in purified 
mtDNA, in X phage titers and the final phage clones 
revealed no statistically significant variation (Chi-square 
test), indicating that the method is unbiased toward which 
molecule is cloned (Table 1). 

Determination of the X phage replication error rate 

We experimentally tested the X phage replication error 
rate by replica-plating a single X clone containing 
mtDNA 60 times in 6 parallel lines, corresponding to 
360 passages of a single A-mtDNA clone. The entire 
mtDNA from the initial clone as well as all final 
six clones were fully sequenced. In comparison with the 
input sequence no new mutations were found, thus 
supporting an extremely low error rate during X phage 
DNA replication, previously estimated to be 7.7 x 10" 8 
mutations/bp (18). 

No evidence for germline recombination of mammalian 
mtDNA 

After having adapted and validated the direct X phage 
method to clone mtDNA from mouse tissue samples, we 



applied this technique to experimentally test whether re- 
combination of mtDNA occurs in the mammalian 
germline. The NZB-BALBc mouse line was generated 
by cytoplast transfer between an NZB-derived cytoplasm 
and a BALBc embryo >20 years ago (14). Heteroplasmy 
levels in these mice segregate in a random manner, best 
explained by random genetic drift and the heteroplasmy 
levels in oocytes are dependent on the levels in the mother 
(14). Different estimates have found that the proportion 
of NZB mtDNA in oocytes ranges from 0.9% to 41% 
(14,19). As these mice have retained an intracellular 
heteroplasmy between NZB and BALB/c mtDNA for 
>20 years, they pose an ideal model to investigate 
mtDNA recombination and clonality in the maternal 
mammalian germline. 

We used the above mentioned RFLP analysis to screen 
single mtDNA molecules cloned in the X phage (Figure 2a 
and b). After screening a total of 2922 mtDNA molecules 
from liver and kidney samples from NZB-BALB/c 
animals (n = 2), we identified only a single aberrant 
molecule with a mixed BALB/c-NZB digestion pattern 
(Figure 2a). Sequencing this entire chimeric molecule 
revealed that the informative SNPs were almost exclu- 
sively NZB-derived with only one single nucleotide 
position representing a BALB/c SNP (15 657c>T). All 
adjacent polymorphisms were NZB (positions 15 603 and 
15 917, Supplementary Table S4), suggesting that this 
variant represents homoplasy, i.e. a single mutational 
event, rather than a recombination event. In total, we 
screened 70128bp of mtDNA by RFLP analysis and 
identified only a single sequence variant, which suggests 
that the mutation rate is ~1.43 x 10" 5 mutations/bp. This 
is in good concordance with the mutation load estimate 
obtained by PCR-clone-sequencing (Supplementary Table 
S5). In order to identify any additional spontaneous or 
homoplastic variants that might have fixed in the two 
genomes, we sequenced the entire mtDNA of nine 
BALBc and eight NZB-derived X clones. All NZB 
mtDNA clones analyzed differed from previously 
reported NZB sequence (20) at two positions (2767g>C, 
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7778g>T), whereas both the NZB and BALB/c molecules 
differed by a single adenine insertion in the polyadenine- 
stretch of the tRNA Arg gene. As we were unable to deter- 
mine whether these changes were present in the initial 
founder animals, we cannot determine whether the NZB 
variants were acquired during the last 20 years of mouse 
breeding and clonally expanded, or whether these variants 
were already present in the original NZB strain. 



DISCUSSION 

The creation of novel recombinant sequence variants 
requires not only an enzymatic process for mtDNA 
recombination, but also the presence of different 
sequence variants of mtDNA. To some extent the exclu- 
sive maternal transmission of mtDNA provides a natural 
barrier against recombination. However, there are reports 
of low-level heteroplasmy in mammals, including humans 
(10,11), suggesting that this may in fact be the prevailing 
situation. Sperm could theoretically provide a source of 
divergent mtDNA and transmission of paternal mtDNA 
occurs in interspecific crosses (12,13). However, under 
normal physiological conditions, paternal mitochondria 
are actively targeted and destroyed in the oocyte by an 
incompletely understood mechanism, which is present 
throughout the animal kingdom (21-24) with a few excep- 
tions (12,13). Recent high-resolution STED microscopy 
studies of primary cells and cell lines in tissue culture 
have revealed that mammalian mitochondrial nucleoids 
typically contain a single copy of mtDNA (25,26), 
further reducing the chance that two molecules will be 
within close proximity to favor recombination events. 
However, the nucleoid organization of mtDNA in the 
maternal germline has not yet been studied and may 
possibly be different from cells in culture. 

Despite arguments against mtDNA recombination, 
there are nevertheless many reports suggesting the occur- 
rence of recombination in the mammalian mtDNA. 
Analysis of recombination by PCR amplification is 
prone to artifacts and these results are therefore contro- 
versial (21). In this study, we sequenced cloned PCR 
products and could unequivocally demonstrate relatively 
high levels of polymerase jumping between mtDNA and 
nuclear pseudogene templates, although the latter were 
present at very low levels. We thus conclude that PCR 
amplification prior to cloning is an unreliable method 
when studying mtDNA recombination. These results 
lead us to establish a direct cloning approach by using 
the k phage to trap individual mtDNA molecules prior 
to amplification. Control experiments showed that this X 
phage cloning approach is unbiased and does not create 
artificial recombinants. 

The NZB-BALB/c mouse line analyzed in this study has 
been heteroplasmic for two essentially neutral mtDNA 
variants, that differ at 90 nucleotide positions, for more 
than 50 generations. We applied a direct cloning approach 
followed by sequencing of individual mtDNA molecules 
and were unable to detect any mtDNA recombination 
events despite the extensive period of heteroplasmy. The 
mtDNA germline bottleneck effect allows only a small 



fraction of the mtDNA copies in the maternal germline 
to contribute to the mature oocyte and sequence variants 
can therefore be rapidly fixed within a population (27). 
One previous study in heteroplasmic mice, using a very 
limited number of generations, found no recombination 
in neutral mtDNA variants, whereas <1% recombinant 
molecules were found in a mouse line carrying deleted 
mtDNA (28). Importantly, it should be noted that the 
analysis of these mouse lines mainly covered somatic 
recombination events, but the results nevertheless 
suggest that recombination may be possible under scen- 
arios of mitochondrial dysfunction. The results in our 
study clearly show that recombination is vanishingly 
rare both in somatic tissues and in the germline under 
conditions of neutral heteroplasmy. As low-level 
heteroplasmy under neutral conditions seems to be the 
common state in nature, recombination should not 
contribute to genetic variation and distribution of 
mtDNA variants in most natural animal populations. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 
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